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'"7^ Abstract. Pairing plays a crucial role in nuclear spectra and attempts to describe it has a long 

history in nuclear physics. The limiting case in which all single particle states are degenerate, 
but with different s-wave pairing strengths was only recently solved. In this strong coupling limit 
the nuclear pairing Hamiltonian also exhibits a supersymmetry. Another solution away from 
those limits, namely two non-degenerate single particle states with different pairing strengths, 
was also given. In this contribution these developments are summarized and difficulties with 

JC^ . possible generalizations to more single particle states and d-wave pairing are discussed. 
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f-^ ■ 1. Introduction 

Pairing plays a very prominent role in nuclear structure physics. The s-wave pairing problem 
can be conveniently formulated using the quasi-spin operators: 

S}=J2(-^ U - m) 4m4-^ S-=J2(-l) iJ - m) aj-ma jm , (1) 



m>0 ra>0 



and 
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% = 2 E («] m a i m + a ) -m a i -™ - l ) = **3 - o%' ( 2 ) 



where Qj = j + ^ is the maximum number of pairs that can occupy the level j and the number 
operator is 

^■ = nE ( a ) m a J m + a] _ m a. } _ m ) . (3) 

m>0 

Quasi-spin operators generate mutually commuting SU(2) algebras: 

[S+, Sr] = 25^1 [Si Sf] = ±5 lJ Sf (4) 

which are realized in the Sj = ^j representations. As we describe later, it is possible to define 
similar operators for d-wave pairing situations. 



2. Exact Solutions 

Over the years considerable attention was paid to exactly solvable pairing Hamiltonians with 
one- and two-body interactions. These cases include 

• The exact quasi-spin limit [lj: 

H = -\G\J2S+S7,. (5) 

jj' 

• Richardson's solution for the case when the single particle energies are added to the 
Hamiltonian in Eq.([5]) [2 

H = J2 e i a ) m a j m - \G\ Y, S+SJ, . (6) 

3™, jj' 

• A model of Gaudin's [3], which is closely related to the Richardson's solution. 

• The limit with separable pairing in which the energy levels are degenerate (the one-body 
term becomes a constant for a given number of pairs) [H [5j|6]: 

H = -\G\J2c*c f S+S-,. (7) 

jj' 

• Most general separable case with two orbitals [7]. 

Let us first examine the degenerate case. It is convenient to define the operators 

S+ (0) = J2 c*j $? and S~ (0) = £) Cj Sj . (8) 

j J 

Then the Hamiltonian in Eq. (|7|) can be written as 

H = -\G\S + (0)§-(0). (9) 

In the 1970's Talmi showed that under certain assumptions, a state of the form 

5+(0)|0)=]T C *5+|0) (10) 

j 

where |0) is the particle vacuum, is an eigenstate of a class of Hamiltonians including the one 
above. Indeed 

HS+(0)\0) = (-|G|£%| C /j 5+(0)|0). (11) 

How can we identify other one-pair states, which are eigenstates of the Hamiltonian in Eq. ((7])? 
For example for two levels j\ and j'2, the state orthogonal to the one in Eq. (fT0|) 

is also an eigenstate with E=0. It turns out that there is a systematic way to find such states. 
First one introduces the operators 

S + W = E , ,l„J t, S-(x) = V— S-ra-S". (13) 
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Cj\ 2 X 3 



(Note that the operators in Eq. ([8]) are the same as the operators in Eq. (fT3j) calculated at 
x = 0). Then one can show that [4, 5j 

S+(0)S+(z[ N) )... £+(4^)10) (14) 

is an eigenstate of the Hamiltonian in Eq. ([7]) with energy 

^ = -|gi(e%i^i 2 -E4)] ( 15 ) 

V j fc=i % / 

if the following Bethe ansatz equations are satisfied: 

E ~? J ' /2 W =-rU+ E IF^T m = l,2,...iV-l. (16) 



1/U.I2 _W _W ' f- W) W) 

J L /\ c j\ z 'm z ™ k=l(k^m) z m z k 



Similarly 



S+(xW)S+( X W)...S+( X W)\0) (17) 

is an eigenstate with zero energy if the following Bethe ansatz equations are satisfied: 

-0 72 N 1 

E 17 , ./_ (TV) = E (N)_ (N) for every m = 1, 2, . . . ,N. (18) 

j Vl c il x *n k=l(k^m) x m x k 



It should be emphasized that the states in Eqs. (|14p and (|17p are eigenstates of the Hamiltonian 
in Eq. (J7|) if available single-particle levels are at most half full. One can show that, if the 
single-particle levels are more than half full, the state 

S-{z^)S^z^)...S-{&m (19) 

is an eigenstate with the same energy as in Eq. (|15p if the Bethe ansatz equations given in 
Eq. (fT6|) are satisfied [5]. In Eq. (fT9|) |0) designates the state where all single-particle levels are 
completely filled. 

One should emphasize that although the solutions presented here are exact solutions, to 
obtain explicit expressions for the energies and eigenstates, one still needs to obtain solutions of 
the Bethe ansatz equations. So far methods to solve the Bethe ansatz equations were developed 
only in a limited number of cases [SJ [9]. 

Note that it is also possible to calculate the quantum invariants of these pairing Hamiltonians 

GDI. 



3. Supersymmetry 

The above discussion indicates that the states of the Hamiltonian in Eq. (|7|) with N pairs (Eq. 
T4|) and with N max + 1 — N pairs (Eq. [T7|) have the same eigenenergy. Furthermore the zero 
energy states are missing if the single-particle states are more than half full. This situation is 
reminiscent of the spectral conditions of the super symmetric quantum mechanics |llj . Indeed 
introducing the operator [6] 



f = exp Lil E(^ + an) 



(20) 



which transforms the empty single-particle state, |0), to the fully occupied state, |0): 

f|0) = |0), (21) 

one can define new operators 

B-=f*§-(0), B + = S + (0)f. (22) 

Supersymmetric quantum mechanics tells us that the partner Hamiltonians Hi = B + B~ and 
H2 = B~B + have identical spectra except for the extremal (usually ground) state of H\. Here 
two Hamiltonians H\ and 62 are actually identical and equal to the pairing Hamiltonian, Eq. 
0. Hence the role of the supersymmetry is to connect the states in Eqs. (j!4j) and (|19p . i.e. this 
supersymmetry connects particle and hole states. 

4. Exact solution for two single-particle states 

It turns out that one can find an exact solution for the case where there are only two single- 
particle levels [TJ: 

ill = E M° - E <$<*%% + E e i%> (23) 

3 33' 3 

where £j and Cj's are dimensionless and the sums are performed over only two single-particle 
states. 

The eigenstates of the Hamiltonian in Eq. (I23p can be written using the step operators: 



^ + (*) = E 9 . °\ , 2 Sj (24) 

*-? 2,£j — \Cj\ z x J 

as 

J + (x 1 )J + (x 2 )...J + (x N )\0). (25) 

Defining the auxiliary quantities 
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/3 = 2 , ^1 ^ ,5 = 2 J2J1 ' , Jll ,f , (26) 

l°Jll l t -J2 I l C Jll 1^2 I 



one obtains the energy eigenvalues as 



N r 

OXr 



n=l P Xn 

In the above equations, the parameters x\~ are to be found by solving the Bethe ansatz equations 

O-lr-l 2 8 N 2 

y^ "jI c jI _ P _|_ y^ L C28) 

j 2 £j — |Cj-| Xfc p — Xk n=1 /_^ fc x ^n — ^fc 

5. Prospects for other exact solutions 

We have shown that there is an exact solution of the s-wave pairing problem with two non- 
degenerate orbitals in terms of the solutions of the Bethe ansatz equations. (This problem is of 
course numerically diagonalizable in an SU(2) x SU(2) basis). It seems to be very difficult to 
generalize the Bethe ansatz method to the case of three non-degenerate orbitals. (Of course the 
three single-particle level problem is diagonalizable in an SU(2) x SU(2) x SU(2) basis). 



However, the Bethe ansatz method seems to be generalizable to at least some d-wave pairing 
situations. In his work, mentioned earlier Talmi showed that, for a Hamiltonian with only one- 
and two-body interactions, if the doubly-magic ground state energy is normalized to zero, 

H\0) = 0, (29) 

then the conditions 

[H,S\0)] = VS\0) (30) 

and 

1[H,S\0)],SH0)} = W(SH0)) 2 (31) 

are sufficient to determine the entire energy spectra. In addition, if one considers d-wave pairing 
operator 

33' 

then the additional condition 

[[H,S%Dl I ] = WS^D\ I (33) 

makes eigenstates with d-wave pairs possible 

HDl\0)<xDl\0). (34) 

Note that these conditions suggest existence of an algebraic structure. In fact, Ginocchio model 
[12j satisfies the double commutators given above. 

One should explore if the Bethe ansatz method can be generalized in other cases where we 
also have an algebraic framework. A possible answer lies in Gaudin's method mentioned earlier. 
It can be shown that the infinite-dimensional algebra 

[J+(A),J- M ] = 2 J0(A ]- J0 ^ , (35) 

A — fi 

[J°(\),J±(»)}=± J±iX) ~ J±M , (36) 

A — [I 

[J°(A),J°( M )] = [J ± (A),J ± ( / z)]=0 (37) 

can be used to find eigenvalues and the eigenvectors of the Hamiltonian 

H(X) = J°(A)J°(A) + V(A)J-(A) + V(A)J+(A). (38) 

For details the reader is referred to, e.g. Refs. [13], [H], and [15]. (It should be emphasized that 
the Hamiltonian in Eq. f)38f) is not the Casimir operator of the infinite-dimensional algebra given 
above). The connection between the Hamiltonian in Eq. (j38|) and the Hamiltonian originally 
solved by Richardson (Eq. [6]) is uncovered by using the realization 

j0 ( A ) = Er^r and J± ( A ) = £r^T' (39) 

1=1 % 1=1 l 

where Sf and Sf are the quasi-spin operators. A possible procedure to include d-wave pairing 
would first identify the underlying algebraic structure (i.e. the SO (8) symmetry of the Ginocchio 
Model), then introduce the analogs of the Eqs. ([35]) through (|38|) using a realization analogous 



to that in Eq. (f39|h In fact a boson pairing model, using the SU(1,1) analog of the quasi-spin in 
this fashion, was already studied in some detail |144ll6j. However, such a program with arbitrary 
fermion pairings remains to be largely unexplored. 

One should finally remark that some of the techniques discussed here are more generally 
applicable to other Hamiltonians with special one- and two-body interactions. One such example 
is the Hamiltonian describing a dense electron gas near the center of a core-collapse supernova 
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